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BACKGROUND OF THE INVENTION 
[0004] Pursuant to the MPEG-2 Advanced Audio Coding (MPEG-2 
AAC) standard, audio signals are sampled frequencies 
starting from 8 IcHz to 96 JcHz . The samples are grouped into 
consecutive frames of 1024 samples. The frames are grouped 
into windows that comprise 2048 samples. However, each 
window has a 50% overlap with the previous window. 
Accordingly, the first 1024 samples of a window are the 
same as the last 1024 samples of the previous window. A 
window function is applied to each window, resulting in 
sets of 2048 windowed samples. The modified discrete cosine 
transformation (MDCT) is applied to each set of windowed 
samples, resulting in 1024 frequency coefficients. The 
frequency coefficients are then quantized and coded for 
transmission. 
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[0005] The first step in decoding is to establish the frame 
synchronization. Once the frame synchronization is found, 
the AAC bitstream is demultiplexed. This includes Huffman 
decoding, scale factor decoding, and the decoding of the 
side information used in tools such as mono/stereo, 
intensity stereo, TNS, and the filter bank. The spectral 
samples are decoded and copied to the output buffer in a 
samples fashion. After Huffman decoding, each coefficient 
must be inverse quant izated by a 4/3 power nonlinearity and 
then scaled by the quantizer step size. Finally, the 
Inverse MDCT (IMDCT) transforms the spectral coefficients 
into time domain. After the IMDCT transform, the output 
samples are windowed, overlapped, and added for generating 
the final pulse code modulate (PCM) samples. 

[0006] As each block is encoded independently, the 
effect of quantization noise is different in different, 
blocks. This results in a large noise at the block 
boundaries when the blocks are decoded. The MDCT is a 
linear orthogonal lapped transform, based on the idea of 
time domain aliasing cancellation (TDAC) . Though the MDCT 
is 50% overlapped, a sequence data after MDCT has the same 
number of coefficients as samples before the transform 
(after overlap-and-add) . This means that a single block of 
IMDCT data does not correspond to the original block on 
which the MDCT was performed. When subsequent blocks of 
inverse transformed data are added using 50% overlap, the 
errors introduced by the transform cancels out the TDAC and 
effectively removes an otherwise easily detectable blocking 
artifact between transform blocks, 

[0007] The IMDCT output is mapped to the full spectrum for 
subsequent use in the window overlap add to implement the 
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window overlap add. For a window sequence of length 2048, 
the memory includes a present frame buffer of 2048 samples, 
previous frame buffer of 1024, and a window coefficient 
buffer of 2048 coefficients. The implementation of the 
window overlap add is computationally intensive. 
[0008] Further limitations and disadvantages of 
conventional and traditional systems will become apparent 
to one of skill in the art through comparison of such 
systems with the invention as set forth in the remainder of 
the present application with reference to the drawings. 
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BRIEF SUMMARY OF THE INVENTION 

[0009] Presented herein is a reduced memory implementation 
technique of filterbank and block switching for real-time 
audio applications. In one embodiment, there is presented a 
method for calculating pulse code modulated samples. The 
pulse code modulated samples are calculated by accessing an 
IMDCT sample from a previous set of IMDCT samples, 
accessing an IMDCT sample from a present set of IMDCT 
samples, calculating a first pulse code modulated sample 
from the accessed IMDCT sample from the previous set of 
IMDCT samples and the accessed IMDCT sample from the 
present set of IMDCT samples, and calculating a second 
pulse code modulated sample from the accessed IMDCT sample 
from the previous set of IMDCT samples and the accessed 
IMDCT sample from the present set of IMDCT samples. 
[0010] In another embodiment, there is presented a system 
for calculating pulse code modulated samples. The system 
includes a first address register, a second address 
register, and an arithmetic logic unit. The first address 
register accesses an IMDCT sample from a previous set of 
IMDCT samples. The second address register accesses an 
IMDCT sample from a present set of IMDCT samples. The 
arithmetic logic unit calculates a first pulse code 
modulated sample from the accessed IMDCT sample from the 
previous set of IMDCT samples and the accessed IMDCT sample 
from the present set of IMDCT samples and calculates a 
second pulse code modulated sample from the accessed IMDCT 
sample from the previous set of IMDCT samples and the 
accessed IMDCT sample from the present set of IMDCT 
samples . 
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[0011] In another embodiment, there is presented a circuit 
for calculating PCM samples. The circuit includes a 
processor and an instruction memory. The processor executes 
a plurality of executable instructions. The instruction 
memory stores the plurality of executable instructions. 
Execution of the executable instructions causes accessing 
an IMDCT sample from a previous set of IMDCT samples from a 
first memory, accessing an IMDCT sample from a present set 
of IMDCT samples from a second memory, calculating a first 
pulse code modulated sample from the accessed IMDCT sample 
from the previous set of IMDCT samples and the accessed 
IMDCT sample from the present set of IMDCT samples, and 
calculating a second pulse code modulated sample from the 
accessed IMDCT sample from the previous set of IMDCT 
samples and the accessed IMDCT sample from the present set 
of IMDCT samples. 

[0012] These and other advantages and novel features of the 
present invention, as well as details of an illustrated 
embodiment thereof, will be more fully understood from the 
following description and drawings. 
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BRIEF DESCRIPTION OF SEVERAL VIEWS OF THE DRAWINGS 
[0013] FIGURE 1 is a block diagram describing the encoding 
of an audio signal; 

[0014] FIGURE 2 is a block diagram of an exemplary audio 
decoder in accordance with an embodiment of the present 

invention; 

[0015] FIGURE 3 is a block diagram describing the decoding 
of an audio signal; 

[0016] FIGURE 4 is a block diagram of a circuit in 
accordance with an embodiment of the present invention; and 
[0017] FIGURE 5 is a flow diagram for calculating the pulse 
code modulated samples in accordance with an embodiment of 
the present invention. 
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DETAILED DESCRIPTION OF THE INVENTION 

[0018] Referring now to FIGURE 1, there is illustrated a 
block diagram describing the encoding of an exemplary audio 
signal A(t) . Pursuant to the MPEG-2 Advanced Audio Coding 

(MPEG~2 AAC) standard, the audio signal is sampled at rates 
starting at 8K samples/sec to 96K samples/sec. The samples 
are grouped into frames Fo...Fn of 1024 samples, e.g., 
Fx(0) . . .Fx(1023) . The frames Fo-.-Fn are grouped into windows 
Wo...Wn that comprise 2048 samples, e.g., Wx (0) . . . Wx (2047) . 
However, each window Wx has a 50% overlap with the previous 
window Wx-i- Accordingly, the first 1024 samples of a window 
Wx are the same as the last 1024 samples of the previous 
window Wx-i. A window function w(t) is applied to each 
window Wo...Wn, resulting in sets wWo...wWn of 2048 windowed 
samples, e.g., wWx(O) . . .wWx(2047) . The modified discrete 
cosine transformation (MDCT) is applied to each set 
wWo- . .wWn of windowed samples wWx(O) . . .wWx(2047) , resulting 
sets MDCTo . . . MDCTn of 1024 frequency coefficients, e.g., 
MDCTx(O) . . .MDCTx(1023) . 

[0019] The sets of frequency coefficients MDCTq ... MDCTn are 
then quantized and coded for transmission, forming what is 
known as an audio elementary stream AES, The AES can be 
multiplexed with other AES. The multiplexed signal, known 
as the Audio Transport Signal (TS) can then be stored 
and/or transported for playback on a playback device. The 
playback device can either be local or remotely located. 
Where the playback device is remotely located, the 
multiplexed signal is transported over a communication 
medium, such as the internet. During playback, the Audio TS 
is demultiplexed, resulting in the constituent AES signals. 
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The constituent AES signals are then decoded, resulting in 
the audio signal. 

[0020] Referring now to FIGURE 2, there is illustrated a 
block diagram describing an exemplary audio decoder 3 00 in 
accordance with an embodiment of the present invention. 
Once the frame synchronization is found, the AAC bit stream 
is demultiplexed by a bitstream demultiplexer 305. The 
bitstream demultiplexer separates the parts of the MPEG-2 
AAC data stream into the parts for each tool, and provides 
each of the tools with the bitstream information related to 
that tool. The AAC decoder includes Huffman decoding 310, 
scale factor decoding 315, and the decoding of the side 
information used in tools such as mono/stereo 32 0, 
intensity stereo 325, TNS 330, and the filter bank 335. The 
sets of frequency coefficients MDCTq. . .MDCTn are decoded and 
copied to the output buffer in a sample fashion. After 
Huffman decoding 310, an inverse quantizer 34 0 inverse 
quantizes each set of frequency coefficients MDCTq. . -MDCTn 
by a 4/3 power nonlinearity . The scale factors 315 are 
then used to scale sets of frequency coefficients 
^ MDCTq ... MDCTn by the quantizer step size. Additionally, 
tools including the mono/stereo 320, intensity stereo 325, 
TNS 33 0, and can apply further functions to the sets of 
frequency coefficients MDCTq. . .MDCTn- Finally, the filter 
bank 335 transforms the frequency coefficients MDCTq ... MDCTn 
into the time domain signal A(t) . 

[0021] The filter bank 335 transforms the frequency 
coefficients by application of the Inverse MDCT (IMDCT) , 
the inverse window function, window overlap, and window 
adding. Referring now to FIGURE 3, there is illustrated a 
block diagram describing the transformation of the 
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frequency coefficients MDCTq . . - MDCTn into the time domain 
signals A{t) . Application of the IMDCT to the sets of 
frequency coefficients MDCTq. . .MDCTn results in overlapping 
sets IMDCTo . . . IMDCTn of 2048 IMDCT samples, e.g., 
IMDCTx(O) . . . IMDCTx (2047) . The inverse window function w"^(t) 
is applied to each set of IMDCTo ... IMDCTn of 2048 IMDCT 
samples IMDCTx (0) ... IMDCTx (2047) , resulting in overlapping 
sets w'^ IMDCTo- . -w"^ IMDCTn of 2048 dewindowed samples, e.g., 
w"^IMDCTx{0) . . .w'^iMDCTx(2047) . The overlapping sets w" 
^IMDCTo .w'^ I MDCTn of 2048 dewindowed samples w" 
^IMDCTx(O) . . .w~^IMDCTx(2047) are then added, e.g., w" 
^IMDCTx(O) . . .w'^IMDCTx(1023) is added to w'^^IMDCTx-i ( 1024 ) • - -w^ 
^IMDCTx-i (2047) , e.g., resulting in the frames F(0)...F(n) 
with 1024 PCM samples Fx (0) . • . Fx (1023) • 

[0022] The symmetry of the IMDCT output can be exploited to 
simplify the operation of applying the inverse window w*'^(t) 
to .each set of IMDCT© ... IMDCTn of 2048 IMDCT samples 
IMDCTx (0) - . - IMDCTx (2047) / and adding the overlapping 
portions of the dewindowed samples, e.g., w'^IMDCTx (0) . . . w" 
^IMDCTx ( 1023 ) , and w"^IMDCTx-i (1024) .w"^IMDCTx-i (2047) , As a 
result, there is no need to unfold the IMDCT output 
IMDCT (0) ... IMDCT (n) to the full spectrum and the samples 
are read only once. The foregoing results in the reduced 
memory requirements and processor operations. 
[0023] For 2048-point IMDCT, it is noted that IMDCTx (0) = - 
IMDCTx(1023), IMDCTx(l) = -IMDCTx (1022) , and IMDCTx(1024) = 
IMDCTx(2047) , and IMDCTx(1025) = IMDCTx (2046) . The 
dewindowed samples w'^IMDCTx (0) . . . w"^IMDCTx(n) can be 
calculated with the following formula: 

w'^IMDCTx(i) = w"^(x) * IMDCTx (i) 

where w'^{i) is a constant 
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[0024] The samples Fx (0) . . . Fx (1023 ) for frame Fx are 
determined by adding the overlapping portions of the 
dewindowed samples w'^IMDCTx (0) . . . w'^IMDCTx (1023 ) , and w" 
^IMDCTx-i (1024) .w"^IMDCTx-i (2047) . As a result, for frame 
Fx, the samples Fx(0) . . .Fx(1023) can be calculated with the 
following formula: 

Fx(i) = w(i)*IMDCTx(i) + 

w(i + 1024) *IMDCTx-i (i + 1024) 

where i = 0 to 1023 

X is the present window index 
x-1 is the previous window index 

[0025] A brute force method would require storage of IMDCTx- 
i(1024) . , .IMDCTx-i(2047) , IMDCTx(O) . . . IMDCTx(1023) , and 

w(0) . - .w(2047) , a total of 4096 memory locations, as well 
as 1024 locations for Fx(0) . . .Fx(1023) . Additionally, two 
memory accesses from the IMDCT samples and two accesses 
from the inverse window function w"^ (x) would be needed for 
each sample Fx (0) . . . Fx (1023) , for a total of 4096 memory 
accesses. 

[0026] The symmetry of the IMDCT samples 

IMDCTx(O) . . .IMDCTx(2047) can be exploited to reduce the 
total number of memory locations and memory accesses needed 
to calculate the samples Fx(0) . . . Fx(1023) from IMDCTx- 
i(1024) . . . IMDCTx-i (2 047) , IMDCTx(O) . . . IMDCTx ( 1023 ) . 
[0027] As noted above, IMDCTx(x) = -IMDCTx (1023 -x) , for x = 

0 to 1023, and IMDCTx- 1 (x) = IMDCTx-i (3071-x) , for x = 1024 
to 2047. As a result, IMDCTx (0) ... IMDCTx (511) , and IMDCTx- 

1 (1024) ... IMDCTx-i (1535) can be stored, while 
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IMDCTx{512) . . .IiyiDCTx(1023) , and IMDCTx-i (1536) . . . IMDCTx- 

i(2047) can be determined from IMDCTx (0) . . . IMDCTx (511) , and 
IMDCTx-i (1024) . . . IMDCTx-i (1535) , respectively. The foregoing 
results in reduced memory consumption. 

[0028] Additionally, the number of memory accesses can also 
be reduced. For i = 512 to 1023, Fx(i) can be calculated by 
the following formula: 

Fx(i) = w(i) *-IMDCTx(1023-i) 
+ w(i + 1024) *IMDCTx-i(2047-i) 
where i = 512 to 1023 

For i = 0 to 511, Fx(i) is calculated as: 

Fx(i) = w(i) *IMDCTx(i) + 

w(i + 1024) *IMDCTx-i (i + 1024) 
where i = 0 to 511 

[0029] As can be seen, Fx(0) and Fx (1023) are both 
calculated from IMDCTx(O) and IMDCTx-i (1024) . Similarly, for 
every i from 0 to 511, Fx(i) and Fx(1023-i) are calculated 
from the same IMDCT values, IiyiDCTx(i) and IMDCTx-i (i4-1024 ) . 
Accordingly, Fx(i) and Fx(1023-i) can be calculated from 
the same memory access from the stored IMDCT samples . 
[0030] The samples Fx(0) . . .Fx(1023) can calculated by 
executing the following operations for i = 0 to 511. 

Fx(i) = w(i) *IMDCTx(i) + 

w(i + 1024) *IMDCTx-i (i + 1024) 
Fx(1023 - i) = w(1023-i) *-IMDCTx(i) 

+ w(2047-i) *IMDCTx-i (i + 1024) 
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[0031] As can be seen, for each iteration, only two memory 
accesses are made from the stored IMDCT samples, and four 
memory accesses from the inverse window function samples 
are made, for a total of 3 0 72 memory accesses. The 
foregoing can be used for the samples for frames FO...Fn, 
thereby reconstructing the signal A(t) . 

[0032] Referring now to FIGURE 4, there is illustrated a 
block diagram describing an exemplary circuit for 
calculating PCM samples Fx(0) . . .Fx(1023) for a frame Fx. The 
circuit can form a portion of the gain control 340. The 
circuit includes a PCM buffer 505, a previous window IMDCT 
buffer 510, a present window IMDCT buffer 515, an inverse 
window buffer 520, and a processor 522. The processor 522 
includes address registers 530, and an arithmetic logic 
unit (ALU) 525, The previous window IMDCT buffer 510 
comprises 512 memory locations 510 (1024) ... 510 (1535) for 
storing IMDCT samples, IMDCTx-i (1024) . . . IMDCTx-i (1535) . The 
present window IMDCT buffer 515 comprises 512 memory 
locations 515 (0 ) . . . 515 (511) for storing IMDCT samples, 
IMDCTx(O) . . . IMDCTx-i (511) . The inverse window buffer 520 
comprises 1024 memory locations 520 (0) ... 520 (1023) for 
storing the inverse window buffer 520 (0) ... 520 (1023) . The 
PCM buffer 505 comprises 1024 locations 505 (0) ... 505 (1023) , 
that are associated with PCM samples Fx (0) . . . Fx (1023 ) for a 
frame Fx. 

[0033] The address registers includes address registers, 
wptrl, wptr2, IMDCTx-iptr, IMDCTxptr, PCMBuf TopPtr , and 
PCMBuf BottomPtr . The address registers wptrl, wptr2, 
initially point to inverse window memory locations 520(0), 
and, 520(1023), , respectively. The inverse window memory 
locations 520(0), 520(1024), 520(1023), and 520(2047) store 
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inverse window samples w"^(0), w"^(1024), w'^(1023), and w' 
^(2047) , respectively. 

[0034] The address register IMDCTx-iptr initially points to 
the previous window IMDCT buffer memory location 510(1024). 
The previous window IMDCT buffer memory location 510(1024) 
stores IMDCTx-i (1024) . The address register IMDCTxPtr 
initially points to the present window IMDCT buffer memory 
location 515(0). The present window IMDCT buffer memory 
location 515(0) stores IMDCTx(O). 

[0035] The PCMBufTopPtr initially points to the PCM Buffer 
memory location 505(0). The PCMBuf BottomPtr initially 
points to the PCM Buffer memory location 505(1023). 
[0036] The contents of the memory locations that are 
referenced by the address registers, wptrl, wptr2, IMDCTx- 
iptr, and IMDCTxPtr are accessed and the ALU 52 5 performs 
arithmetic operations thereon. The results of the 

arithmetic operations are written to the PCM buffer 505 
memory locations referenced by PCMBufTopPtr and 
PCMBuf BottomPtr . 

[0037] The result of the following arithmetic operation is 
stored at the memory location reference by PCMBufTopPtr: 

PCMBufTopPtr = wptrl*IMDCTxPtr + wptr2*IMDCTx-iptr 

[0038] The result of the following arithmetic operation is 
stored at the memory location referenced by 
PCMBuf BottomPtr : 

PCMBufBottomPtr = -wptr2*IMDCTxptr + wptrl*IMDCTx-iptr 
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[0039] After the following operations are performed, 
PCMBufTopPtr, wptrl, IMDCTxPtr, and IMDCTx-iptr are 
incremented, while PCMBuf BottomPtr and wptr2 are 

decremented. After repeating the foregoing for 512 
iterations, the PCM buffer 505 stores the PCM samples for a 
frame Fx- 

[0040] The pseudo-code for the iterations can be described 
as : 

for (i=0; i<512;, i++) 

{ 

*PCMBufTopPtr++ = *wptrl* (*IMDCTxptr) + 

*wptr2* (*IMDCTx-iptr) 

*PCMBuf BottomPtr-- = - (*wptr2 - - ) * ( *IMDCTxPtr++) 

+ (*wptrl + + ) * (*IMDCTx.iptr++) 

} 

[0041] Tables 1 and 2 describe the operations that are 
performed for both a brute force implementation and an 
impementation taking advantage of the IMDCT symmetry, 
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[0042] Referring now to FIGURE 5, there is illustrated a 
flow diagram for calculating PCM samples Fx(0) . . .Fx(1023) . 
At 602, the IMDCT samples IMDCTx (0) . . . IMDCTx (1023) from the 
present set IMDCTx of samples are stored in the present 
IMDCT buffer 515. At 605, IMDCT samples IMDCTx-i 
(1024) . . . IMDCTx-i (2047) from the previous set IMDCTx-i of 
IMDCT samples are stored in the previous IMDCT buffer 510. 
[0043] At 607, the IMDCTx-iptr is set to point at the 
previous IMDCT buffer location 510(1024) storing IMDCTx-i 
(1024) . At 610, the IMDCTxPtr is set to point at the present 
IMDCT buffer location 515(0) storing IMDCTx (0) . 
[0044] From 612-615, the inverse window address 
registers wptrl, wptr2, are initialized. At 612, wptrl is 
set to point to the inverse window buffer location 520(0) 
storing inverse window coefficient w"^(0). At 615, wptr2 is 
set to point to the inverse window buffer location 
520 (1023) storing inverse window coefficient w""''(1023). 
[0045] At 622-625, PCMBufTopPtr and PCMBuf BottomPtr are 
initialized. At 622, PCMBufTopPtr is set to point to the 
PCM buffer locat ion 505(0) that is associated with Fx(0) . At 
625, PCMBuf BottomPtr is set to point to the PCM buffer 
location 505(1023) that is associated with Fx(1023). 
[0046] At 630-634, the contents at the locations 
referenced by the address registers IMDCTx-iptr (630) , 
IMDCTxptr (632), wptrl, wptr2 (634), are fetched. 
[0047] At 640-647, the first PCM value is calculated 
from the fetched contents from the locations referenced by 
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the address registers IMDCTx-iptr , IMDCTxptr, wptrl, wptr2 . 
At 640, the contents of the location referenced by IMDCTx- 
iptr are multiplied by contents of the location referenced 
by wptr2 . At 642, the contents of the location referenced 
by IMDCTxptr are multiplied by the contents of the location 
referenced by wptrl . At 645, the product during 640 during 
642 are added. The foregoing sum is the first PCM value. At 
647, the first PCM value is written to the memory, location 
in the PCM buffer 505 referenced by PCMBuf TopPtr . 
[0048] At 650-657, the second PCM value is calculated 
from the fetched contents from the locations referenced by 
the address registers IMDCTx-iptr, IMDCTxptr, wptrl, wptr2 . 
At 650, the contents of the location referenced by IMDCTx- 
iptr are multiplied by contents of the location referenced 
by wptrl. At 652, the contents of the location referenced 
by IMDCTxptr are multiplied by the contents of the location 
referenced by wptr2 . At 655, the product during 650 is 
subtracted from the product during 652 . The foregoing 
difference is the second PCM value. At 657, the second PCM 
value is written to the memory location in the PCM buffer 
505 referenced by PCMBuf BottomPtr . 

[0049] At 660, the address registers PCMBuf TopPtr , 
wptrl, IMDCTxptr, and IMDCTx-iptr are incremented, while at 
665, the address registers PCMBuf BottomPtr, wptr2 are 
decremented. At 670, a determination is made whether 
PCMBufTopPtr points to location 505(512). If during 670 
PCMBufTopPtr points to location 505(512), the calculation 
of the PCM samples Fx ( 0 ) . . . Fx ( 1023 ) is completed for Fx and 
602-670 are repeated for the PCM samples Fx+i (0) . . . Fx+i (1023) 
for frame Fx+i. If during 670, PCMBufTopPtr does not point 
to location 505(512), then 630-670 are repeated. 
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[0050] The systems and circuits as described herein may 
be implemented as a board level product, as a single chip, 
application specific integrated circuit (ASIC) , or with 
varying levels of the decoder system integrated with other 
portions of the system as separate components. The degree 
of integration of the decoder system will primarily be 
determined by speed and cost considerations. Because of 
the sophisticated nature of modern processor, it is 
possible to utilize a commercially available processor, 
which may be implemented external to an ASIC 
implementation. Alternatively, if the processor is 
available as an ASIC core or logic block, then the 
commercially available processor can be implemented as part 
of an ASIC device wherein certain operations are 
implemented in firmware. 

[0051] While the invention has been described with 
reference to certain embodiments, it will be understood by 
those skilled in the art that various changes may be made 
and equivalents may be substituted without departing from 
the scope of the invention. In addition, many modifications 
may be made to adapt particular situation or material to 
the teachings of the invention without departing from its 
scope. Therefore, it is intended that the invention not be 
limited to the particular embodiment (s) disclosed, but that 
the invention will include all embodiments falling within 
the scope of the appended claims. 
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